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The continuum limit of the tilted SU(2) spin model is shown to give rise to the gauge Landau- 
Lifshitz equation which provides a unified description for various spin orders. For a definite gauge, 
we find a double periodic solution, where the conical spiral, in-plane spiral, helical, and ferromagnetic 
spin orders become special cases, respectively. For another gauge, we obtain the skyrmion-crystal 
solution. By simulating the influence of magnetic field and temperature for our covariant model, we 
find a spontaneous formation of skyrmion- fragment lattice and obtain a wider range of skyrmion- 
crystal phase in comparison to the conventional Dzyaloshinsky-Moriya model. 
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There has been spectacular progress in the study on 
the magnetoelectric effects, which is expected to make a 
realistic step toward an electrical control of magnet ism [Ti- 
ll]. Within the intertwining of theory and experiment, a 
mechanism based on spin-current showed that the en- 
hanced ferroelectric domains can be realized through cy- 
cloidal and conical spin states in certain materials [2j[5]. 
For example, spiral spin state was shown [5 to cause 
electronic polarization. Moreover, complex spin texture 
is interesting by its own right. Recently, skyrmion lattice 
is observed in bulk MnSi [6] and thin film Fe^Coi-^Si [7]. 
Magnons in helical magnets are also probed by neutron 
scattering experiments [8]. Those nontrivial spin tex- 
tures, undoubtedly, play an important role in novel mul- 
tiferroic materials. There are mainly two kinds of mech- 
anisms that cause spiral spin state, one is ferro/antiferro- 
magnetic exchange competition that is believed to be 
the origin of spiral state in manganites [9j, the other is 
the antisymmetric Dzyaloshinsky-Moriya (DM) interac- 
tion [TQl [TT] which arises from spin-orbit (SO) interac- 
tion and manifests in crystal without inversion symme- 
try like MnSi. Historically, Moriya [11] was the first to 
give a microscopic treatment of DM interaction based on 
Anderson's superexchange mechanism with SO interac- 
tion. Thirty two years later, Shekhtman [12] found that 
Moriya's theory has a bond-isotropic form if only one 
takes all the terms up to second order in SO interaction 
which is usually neglected until present. It is important 
to set up a unified description for various spin orders. 

In this letter, we indicate that such a system can be 
described by tilted Heisenberg model in which the tilting 
is related to the effects associated with bonds. We for- 
mulate gauge Landau-Lifshitz equation from this model 
and find solutions of various spin textures and derive the 
dispersion relation of the relevant spin waves, which pro- 
vides a unified theory for spin orders with insight in the 
gauge and geometric point of view. Then we investigate 
the influence of external magnetic field and temperature 
and plot the corresponding phase diagram by making use 



of Monte Carlo simulations. 

In order to reach a unified description of various spin 
ordered phases including the situations beyond the tradi- 
tional ferromagnetic one, we consider a much more gen- 
eralized Heisenberg Hamiltonian 

H = -J J] UjSpriUySpUy 1 (1) 

(33 f ),c 

where c = 1,2,3, j G L with L the lattice space, and 
Sj denotes the c-th component of spin operator at site 
j. These spin operators, proportional to the infinites- 
imal generators of SU(2), obey [S%, Sf] = ihSjie abc S^ 
that governs the time developments of any observable via 
Heisenberg equation of motion for a definite model ([!]). 
In Eq. ([!]), (jf) means the summation is taken over the 
nearest neighbor lattice sites, and the local tilting field Uj 
accounts for any effects arising from either (both) compli- 
cated crystalline fields or (and) cumbersome charge order 
in whatever intricate materials. 

As there exists a homomorphism between SU(2) and 
SO (3) Lie groups, UjSjU' 1 = SjOj in which Sj de- 
notes (Sj,Sj,Sj) and Oj the representation of SO(3), 
each nearest-neighbor term in Eq. ([I]) can be rearranged, 

i.e., Y.<'; s j<\, l <'r s r r r : = SjOj(S r o r ) T = 

c 

SjOjOj, 1 Sj, . Here Oj, = 0~} for orthogonal group has 
been used. Because j' is close to j when the lattice con- 
stant a is taken as an infinitesimal parameter, we can 
expend OjOj, 1 in the vicinity of identity, namely 

OjOJ, 1 = l-aAt(j)i c , (2) 

where we have considered the coordinate of site j' is sim- 
ply that of j plus a bond vector ae v in which e v refers to 
the unit vectors connecting neighborhood of a given lat- 
tice structure. Here £ c denote the representation matrices 
of the infinitesimal generators of SO (3) Lie group, they 
are 3x3 matrices (£ c )ab = tabc and fulfil the commuta- 
tion relations [13 [£ a , iy = —e a b c £ c . Clearly, the feature 
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of the local tilting can be characterized by the SO (3) 
non-Abelian gauge potential = A c v (j)t c which is 

a matrix valued vector field. In order to avoid any am- 
biguity, here we clarify that the j represents a point in 
the lattice space corresponding to the coordinate of real 
space in continuum model, the c labels the component of 
a vector in Lie algebra space while the v labels the one 
in real space. Also for symbol neatness, in Eq. Q and 
thereafter, we write the lattice-site label j of A in paren- 
theses rather than conventional subscripts. By making 
use of Eq. ^ , we can write Eq. ([!]) as 
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Actually, the Casimir invariants Cj = Sj ■ Sj = Sj(sj + 
I) ft 2 in a general system may differ at different lattice 
site, which means the module of spin does not necessarily 
take the same value everywhere. However, in this paper, 
we focus on uniform spin module S in every sites. 

Now we are in the position to make continuum limit, 
—> (l/a) d j d d x, which can be realized by allowing the 
volume per lattice site a d tend to zero and considering the 
lattice label j as a continuous variable r and hence Sj as 
M(r). Equation Q gives rise to the effect Hamiltonian, 



H = 



J 



2a d ~ 



I 



(d u + A v (r)x)M(r) 
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where the additional constant term is omitted. Then 
the corresponding Lagrangian density is given by L = 



- d \M\(cosO-l) 



ji-d 



J/2 (DM) 2 in which! 



refer 



to the azimuthal angles of M. The equation of motion 
for the spin field M(r,i) is derived as the following gauge 
Landau-Lifshitz equation, 



d_ 

dt 



M = a 2 JM x D 2 M, 



(5) 



where D 2 = D V D V and the covariant derivative is given 
by D V M = (d u + A u x)M. Equation ^ is covariant 
under a gauge transformation A^ — >• GA^G -1 +c^GG _1 , 
(M) a -+ E fc G ab (M) b with G g SO(3). 

We first consider a typical gauge field in x-y plane 
A x = (0,0, -gi), A y = (q 2 smq 1 x,-q 2 cosq 1 x,0). We 
find a double periodic solution, Md p , as a steady solu- 
tion of the gauge Landau-Lifshitz equation ([5|, 



' m^x.y) = sm(q 2 y- 
m 2 (x,y) = sin(q 2 y- 
^ m 3 (x,y) = cos(q 2 y- 



■ 0)co8(qix), 
■0) sin(gix), 



(6) 



Here m = (mi, m 2 , ms) refers to M& v /S. This spin order 
is the exact ground-state solution of the system because 
D x Md p = DyM^p = so that the positive definite en- 
ergy functional ^ reaches zero then. Clearly, the conical 



spiral spin order [14] is the special case of q 2 = whose 
special case of /3 = tt/2 reduces to the in-plane spiral spin 
order [15 . The other cases q\ — or q\ = q 2 = corre- 
sponds to a helical spin order [16 or the ferromagnetic 
spin order, respectively. 

To study the excitations above the aforementioned 
ground state (J6|, we take Md p + SM and obtain the fol- 
lowing linearized equation (d t — a 2 JM^ p x D 2 ) SM = 0. 
Since the constraint \M\ = S requires M^ p ■ SM = 0, 
we can assume SM = u(x,y,t)eo + v(x,y,i)e ( / ) with 
the local frame = e z x M^ p /\e z x Md p |, eg = 
x Md p /|e0 x Md p |. Then the equations that possible 
low- lying excitation modes obey are d t u + a 2 JSV 2 v = 
and dtv — a 2 JSV 2 u = 0. Their Fourier transform gives 
rise to the dispersion relation uj 2 = a 4 J 2 S 2 \k\ 4 . One can 
see that the dispersion relation here is happened to be 
the same as that of the spin wave above a ferromagnetic 
ground state in classical Heisenberg model. 

Because the strength tensor F^ y = d x A c y — d y A c x + 
e abc A^A y vanishes for the gauge potential relevant to 
the solution ([6|, the gauge potential can be repre- 
sented as a pure gauge A^ = —G~ 1 d u G with G = 
exp(q 2 y£ 2 ) exp(qix£s). The generating matrix G implies 
an important physical significance, which transforms the 
double periodic spiral order (|6| to the traditional ferro- 
magnetic order, i.e., (Mfe)a = ^Z b G ab (M dp ) b . Here M fe 
is the ground-state solution of Eq. ([5| with null gauge po- 
tential. The double periodic spiral order can be consid- 
ered as a result of parallel displacement of spin with the 
aforementioned gauge potential as connection. Since the 
solutions referring to both orders are in the same equiva- 
lent class of gauge Landau-Lifishitz equation, there would 
be no surprise that the dispersion relations for the exci- 
tations above them are the same. 

Next, we investigate the case with non- vanishing 
strength tensor, which gives rise to skyrmion [17 crys- 
tal solutions. For A x = (-7/J, 0,0), A y = (0,-7/J, 0) 
where 7 denotes the strength of spin-orbit interac- 
tion, we have F xy = (0,0,7 2 /J 2 ) and energy den- 
sity functional: (J/2)d v M ■ d v M + 7M • (V x M) + 
(7 2 /2J) [M 2 + (M 3 ) 2 ]. Here the last term contributes 
an easy-plane anisotropy that is the continuum version 
of Moriya's anisotropic exchange [11], the first two terms 
are the conventional ferromagnetic exchange and the DM 
interaction which was used to explore possible states of 
skyrmion crystal [II [H]. Unlike the solution of double pe- 
riodic spiral order which can be generated through a par- 
allel displacement, we need to solve the gauge Landau- 
Lifshitz equation at present. 

For steady solution of the gauge Landau-Lifshitz equa- 
tion ([5]), it is sufficient to solve the eigen-equation 
D 2 M = XM. The A can be a scalar function in gen- 
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eral while it is assumed to be a constant here for simplic- 
ity. It can be proven that the A is proportional to the 
energy density. As being interested in periodic steady so- 
lution, we can assume M(r) = M (k) ex.-p(ik ■ r). Then 
the eigen-equations become a set of algebraic equations 
for us to determine M(k). We obtain three solutions: 
sinusoidal order with spin paralleling to the wave vector, 
elliptically distorted right-handed and left-handed helical 
order with spin perpendicular to the wave vector. As the 
three eigenvalues depend on | k | 2 merely, we can make su- 
perposition of the eigenmodes corresponding to the same 
eigenvalue. In the present case, 7 > is assumed, so 
we want chose the right-handed helical to construct the 
ground state. The closest-packed lattice of skyrmions 
are superposition of three such eigenmodes with three 
wave vectors of the same length and mutually in 120° 
angle, namely M c = Y^h=i M(ki)e lhi ' r in which k\ = 
(^jfe,|jfe), k 2 = (-^fc,§fc), and k 3 = (0,-fc). Since 
the real and imaginary parts of M c all satisfy Eq. ([5| , we 
can normalize the real part to reach a physical state that 
is a compromise of reducing energy and satisfying the 
unit-length constraint, m s k = Re(M c )/\Re(M c )\. Here 
the unfixed parameter k in m s k determines the lattice 
constant of the skyrmion crystal. 

The particular k is determined by minimizing the aver- 
age energy density which is calculated through numerical 
integration. Some features of the solution is plotted in 
Fig. [T] where spins between the center of skyrmions tend 
to point up although the spins in each skyrmion tend 
to point down. The average energy density of the opti- 
mized configuration of skyrmion crystal is 0.276S 2/ y 2 / J 
with k = O.877/ J, which is higher than helical order's 
0.255' 2 7 2 /J, but the average z-component of spin for 
the solution m s k is +0.17. The skyrmion crystal will 
have lower energy when a sufficient large perpendicu- 
lar magnetic field is applied downwards. Whereas, when 
the magnetic field is further enhanced, a ferromagnetic 
state with the z-component of spin being 1 eventually 
becomes the ground state. This argument is consistent 
with Ref. 0. 

Furthermore, we study what will happen if there exists 
a magnetic anisotropy in the system. Such an anisotropy 
can be introduced by adding the term J2j v(Sj) 2 m the 
spin model in which either the easy axis is chosen as 
z-axis or the easy plane as x-y plane for 77 < or 77 > 0, 
respectively. Choosing this kind of anisotropy is due to 
keeping the original rotational symmetry about the z- 
axis. Then the above gauge Landau-Lifshitz equation 
([5| turns to the following anisotropic one, 
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M = a 2 JM x D 2 M - 2rjM x M', 
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FIG. 1: (Color online) Schematic illustration for spin order of 
skyrmion-crystal solution (left panel) and the corresponding 
distribution of energy density in zero magnetic field (right 
panel). The energy density is lower at the center of skyrmion 
because the local spin chirality is large there that is favored 
by DM interaction; it is larger in the boundaries of skyrmions 
but the spins there favor a perpendicular magnetic field. 



an easy plane anisotropy. When the anisotropy co- 
exist with the aforementioned gauge potential relevant 
to the skyrmion-crystal solution, the eigenequation for 
the original gauge Landau-Lifshitz equation in /c-space 
is modified. One can choose the right-handed helical 
mode (k y , —k x , ip \k\) in which p = y^ 2 + 1 — £, where 
£ = (7 2 + 2Jr]/a 2 ) I (4J7 |fc|). In real space, this mode 
is a helical order of elliptic contour with p referring to 
ratio of semiminor and semimajor axes. It can be seen 
that the larger the 77 is, the smaller the p will be, which 
is in consistent with the requirement for minimizing the 
energy. When £ = we have p = 1, it occurs a cancela- 
tion between the added anisotropy term and the second 
order term of 7 arising from the gauge potential. 

Now we turn to investigate finite temperature effects of 
our covariant model which contains DM interaction and 
magnetic anisotropy simultaneously. For convenience in 
numerical simulation, we start from the lattice version, 



H 



E 



Kb • (S r x S r + ae ) 



(8) 



with M' = (0,0, Ms). Note that the gauge poten- 
tial with skyrmion-crystal solution merely contributes 



where S denotes classical spin of unit module; e refers 
to x or y and r runs through the lattice site of the base 
space; J, K and B denote the exchange, the strength 
of DM interaction and the external magnetic field, re- 
spectively. In our numerical calculation, the Boltz- 
mann constant k B and the lattice spacing a is taken 
as unit. We do Monte Carlo simulations in various 
regimes of model parameters. For zero magnetic field 
B = 0, we find, for a specific strength of DM interaction 
K/J = y/2 tan (27r/6), that the system goes from dis- 
ordered phase to helical phase and then to a new phase 
when temperature is lowering. The new phase (see Fig.[2| 
presents a square lattice of alternatively placed skyrmion 
fragments, some of which appear to be imbedded among 
spin helical textures that is marked by black dot-lines 
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FIG. 2: (Color online) (a) Spin order of sky rmion- fragment 
phase. Here is a patch of spin textures from our simulation 
for K/J = V2tan(27r/6). Along the black dot-line arrows 
the spins form a helical order where each two successive spins 
have relative 60° difference in angle, (b) and (c) The spin z-z 
correlation in the 48 by 48 lattice for various phases in zero 
magnetic field. Temperature is lowering, it goes from disor- 
dered to helical (b) then to sky rmion- fragment (c) phases. 
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FIG. 3: (Color online) Phase diagram of various spin orders 
(Sk means skyrmion lattice, H helical, F ferromagnetic, and 
+ coexistence) in the plane of magnetic field versus temper- 
ature calculated in the 36 by 36 lattice for DM model (left 
panel) and our covariant model (right panel). The color indi- 
cates the total number of skyrmions. The helical phase and 
ferromagnetic phase have no skyrmions and the skyrmion lat- 
tice phase have many skyrmions. The area ratio of skyrmion 
phase to helical phase is larger in our covariant model. 



in Fig. [2] Since the new phase appears in zero mag- 
netic field, it is an emergence of spontaneous formation 
of skyrmion- fragment lattice. For weaker strength of DM 
interaction, saying K/J = y/2 tan (27r/9), we plot the 
phase diagrams in the plane of temperature versus mag- 
netic field based on our Monte Carlo simulations for both 
our model ([8| and the conventional model [7] . Our results 
manifest that the landscape of those two phase diagrams 
are similar while the area ratio of skyrmion lattice phase 
to helical phase in our model is larger than that in the 
conventional model (see Fig. [3| . 

In conclusion, the gauge Landau-Lifshitz equation, as 
the continuum limit of the tilted SU(2) spin model, pro- 
vides a unified description for various spin orders. The 
double periodic solution we found implies the conical spi- 
ral, in-plane spiral, helical, and ferromagnetic spin orders 
as special cases, respectively. The skyrmion-crystal order 
is a solution corresponding to a SO (3) gauge with non- 
vanishing strength tensor. As to the finite temperature 



behavior, a spontaneous formation of skyrmion- fragment 
lattice occurs in zero magnetic field, and the area ratio 
of skyrmion phase to helical phase is larger in our co- 
variant model than in the conventional DM model. Note 
that the magnon band structure observed in recent ex- 
periments [8] does not contradict to double periodic dy- 
namics since it happens when the system is described by 
a three-dimensional gauge potential A x = (—7/ J, 0,0), 
A y = (0,-7/J,0) and A z = (0,0, -7/ J). This gauge 
potential is rotationally invariant and does not bring in 
anisotropy, so our model is equivalent to the conventional 
model and the ground state is the circularly helical state 
with wave vector 7/ J. 
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